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a case study on a 2D triangular lattice in the repulsive Hubbard model 

Hiroaki Ikeda^ *, Yunori NiSiKAWA^ and Kosaku Yamada^ 

^Department of Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502 
"^Synchrotron Radiation Research Center, Japan Atomic Energy Research Institute, Mikazuki, Sayo, Hyogo 679-5148 

(Received February 2, 2008) 

Stimulated by the recent finding of Nao.35Co02.1.3H20 superconductor, we investigate su- 
perconducting instabilities on a 2D triangular lattice in the repulsive Hubbard model. Using 
the third-order perturbation expansion with respect to the on-site repulsion U, we evaluate 
the linearized Dyson-Gor'kov equation. We find that an /-wave spin-triplet pairing is the most 
stable in a wide range of the next nearest neighbor hopping integral t' and an electron number 
density n. The introduction of t' is crucial to adjust the van Hove singularities to the neigh- 
borhood of the Fermi surface crossing around K point. In this case, the bare spin susceptibility 
shows the broad peak around F point. These conditions stabilize the /-wave pairing. Although 
the /-wave pairing is also given by the fluctuation-exchange approximation, the transition tem- 
perature is too low to be observed. This is because the depairing effect by the spin fluctuation 
is over-estimated. Thus, the third-order vertex corrections are important for the spin-triplet 
superconductivity, like the case in Sr2Ru04. 
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Recently, superconductivity in Nao.35Co02.1.3H20 
[ with the transition temperature Tc ~ 5K has been dis- 

■ covered^' and intensively investigated. This compound 
I has a hexagonal layered structure, and the dominant con- 

■ ductive plane is formed by the network of the C0O2 oc- 
tahedra. It is considered as the 2D layered triangular 
lattice. The deintercalation of Na"*" ions corresponds to 

; 0.65/Co hole doping for Co^+{3df with fulfilled isg or- 
. bits in the C0O2 layer. The uniform susceptibility shows 
] the weak-ferromagnetic like behavior, and indicates rel- 
• evance of electron correlations. Up to now, although the 
NMR/NQR measurements has been performed in three 
I groups, they have contradicted each other about the 
pairing symmetry of the superconductivity. However, ab- 
. sence of the coherence peak just below Tc observed by 
Fujimoto et o/.^' has suggested possibility of unconven- 
. tional superconductivity with line-node gap. Quite re- 
' cently, this experimental fact has been confirmed by 
Ishida et al.^^ Thus, we are naturally interested to know 
' what kind of superconductivity can be caused by the cor- 
. related electrons on the 2D triangular lattice. 

About such superconducting instabilities on the tri- 
angular lattice, several theoretical studies have been al- 
ready reported. In RVB superconductivity in the frame- 
work of the t-J model, and in a study using a one- 
loop renormalization group, a d-fid-wave spin-singlet 
pairing has been discussed actively. Possibility of a p-\-ip- 
wave spin-triplet pairing has been proposed phenomeno- 
logically.^^-* Within the paramagnon exchange, a d-wave 
singlet and an odd-frequency s-wave triplet pairings has 
been discussed. In a fluctuation-exchange (FLEX) ap- 
proximation, possibility of an /-wave triplet pairing at 
low carrier limit has been also suggested. 

So far, we have investigated unconventional supercon- 



* E- maii : hiroaki@scphy s . kyot o- u .ac.jp 



ductivity in repulsive Hubbard models for many super- 
conductors including the high-Tc cuprates and Sr2Ru04 
on the basis of the third-order perturbation theory. 
All of them are in agreement with experimental facts, 
especially, concerning the pairing symmetry. Thus, here, 
in the same framework, we investigate superconducting 
instabilities in a simple one-band Hubbard model on the 
2D triangular lattice. In this one-band model, 0.65/Co 
hole doping in Nao.35Co02.1.3H20 corresponds to 1.35 
electron filling, or 0.65 electron filling by the particle- 
hole transformation if the band structure is upside-down. 
Thus, the results for the electron number density n =1.35 
or 0.65 in this model will be helpful for the study of the 
C0O2 superconductor. 

Now, let us consider the repulsive Hubbard model on 
the 2D triangular lattice, 

H = ^^kcl^Cka + U'^n,^n^i. (1) 

k i 

The dispersion relation is given by 

ffc = — 2t I 2 cos cos -t- cos Ky I 

^ ^ (2) 

— 2t y cos V3kx -\- 2 cos — - — cos J ~ A^j 

where t and t' are the nearest neighbor and the next near- 
est neighbor hopping integrals, respectively. The chemi- 
cal potential /i is determined by fixing n. In this system, 
we apply the third-order perturbation expansion with 
respect to U, and evaluate the linearized Dyson-Gor'kov 
(gap) equation. The normal self-energy is given by 

S{k)^J2^^Xo{k~k')go{k') (3a) 

k' 

+ {xlik - k') + ct>l{k + k')) Goik'), (3b) 
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Fig. 1. The contour plot of eigen values for 0.4 < n < 1.0 and 
2 < IZ/t < 8 at t' = and T/t = 0.01. The gray thick line with 
points denotes 1/xo^^'' a-s the reference value of U. The d-wave 
for 0.7 < n < 1.0 and the p-wave for 0.4 < n < 0.6 are stable. 



Fig. 2. The phase diagram of n and t' at U /t = 1.8 and T/t = 
0.01. The eigen values for the black region is above 0.4. The 
/-wave triplet pairing (Si ) is the most stable in the wide region. 



where Xo(g'), <^o(9) = - Qo{p)Qa{ri±p). Since the first- 
order normal self-energy is constant and can be included 
by the chemical potential /x, here we omit it. The gap 
equation is given by 

\A{k) = -J2 U\o{k - k')H-k') (4a) 

k' 

- E ^'xo(fc - k'f {^k') + H-k')) (4b) 

k' 

- E [/32Re [x{q)Qo{k + q)go{k' + q)] J-(fc') (4c) 

k'q 

+ E C/32Re [cl>{q)go{k + q)go{k' + q)] T{-k'), (4d) 

k'q 

for both the spin-triplet and singlet pairs, ^''^ where 
^{k) = \Q{k)\'^ A{k), and A denotes an eigen value for 
each eigen vector A{k). The temperature with A = 1 cor- 
responds to Tc- By estimating the largest positive eigen 
value Amax in the gap equation, we determine which kind 
of symmetry is stable among irreducible representations 
(IR) of the space group D^h: 



IR (symmetry) 


One of Basis funcitons 


B\ (f-wave) 
B2 (/-wave) 
El (p-wave) 
E2 (d-wave) 


sin(^)(cos(^) -cos(^)) 
sin(^)(cos(^) - cos(5^)) 
sin( ) cos( ) 
sin(^)sin(^Ji) 



Less than half — First of all, let us consider the case 
of n < 1.0 less than half filling. In this case, the Fermi 
surface is round, and the density of states is low. Since 
the t'-dependence also is weak, we set t' = 0. In Fig.l, 

we illustrate the contour plot of Amax at T/t = 0.01. The 
white and black regions correspond to Amax < 0.0 and 
Amax > 0.4, respectively. The gray thick line with points 
represents l/x™'^^ as a reference value of U, where Xo^^^ 
is the maximum value of the bare susceptibility Xo{q, 0). 
Thus, actually, U/t may be too large to realize the faded 
region on the right-hand side of this line. Near the half 
filling, the d-wave singlet pairing (£^2) is stable. Far from 
the half filling, it becomes unstable, and possitive eigen 
values are obtained only for small U. Rather, the p-wave 
triplet pairing (£^1) is stable for large U. Thus, in this 



region, the convergence of solutions with respect to U 
is not good. The second-order (Eq.(4a)) and the third- 
order (Eq.(4b-d)) scattering processes prefer different 
pairing states; the former the singlet, the later the triplet. 
However, if assuming that U ~ l/x™'^'^) the d-wave for 
0-''' ^ ^ 1-0 and the p-wave for 0.4 < n < 0.6 are sta- 
ble, and 0.6 < n < 0.7 is on the border line. The electron 
filling n = 0.65 corresponding to Nao.35Co02.1.3Il20 is 
on this border line. Thus, we can not expect the stable 
superconductivity at Tc ~ 5K within this analysis. 

More than half — Next, we consider the case of n > 1.0 
more than half filling. This situation is sensitive to in- 
troduction of t'. Fig. 2 shows the contour plot of Amax at 
U/t = 1.8 and T/t = 0.01. The solid line with Amax = 0.0 
distinguishes four regions. Two d-wave singlet {E2) re- 
gions locate near the half filling n = 1.0. These regions 
are continuously connected with the region of n less than 
half filling. In this case, Amax is less than 0.15 at most. 
On the other hand, far from the half filling, the /-wave 
triplet pairing (Bi) has relatively large Amax, and cov- 
ers a wide region of the phase diagram. Successively, the 
region of p-wave triplet (Ei) spreads only near the low 
carrier n ~ 1.9. This indicates that the /-wave triplet su- 
perconductivity can be realized rather than the d-wave 
singlet one at the carrier number far from the half filling 
on the 2D triangular lattice. This is also the case on the 
density n = 1.35 corresponding to Nao.35Co02.1.3H20. 
The pairing state with Bi symmetry is the most stable 
with large Amax in 0.12 < t' /t < 0.32. Although the pair- 
ing state with E2 symmetry possesses the positive eigen 
value in t' /t > 0.32, we can not expect so high T^. Such 
a difference originates from the location of the van Hove 
singularities near the Fermi surface. Let us discuss the 
situation in details below. 

In Fig.3, we illustrate the dispersion relation and the 
density of states p(uj) at n = 1.35. At t' = 0, the van Hove 
points exist at M point, and K point is the maximum of 
^k- This situation is sensitive to the introduction of t'. 
With increasing t', the energy at K point decreases, and 
that at M point increases. Correspondingly, the van Hove 
points move around K point, and the singularities cross 
the Fermi level at t' /t ^ 0.21. We can also verify it by t'- 
dependence of the peak position in p{uj). Since the Fermi 
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Fig. 3. (a) The dispersion relation ^j^. along the symmetry line 
and (b) The density of states p{w) at n = 1.35 for t' /t = 0.00, 
0.10, 0.15, 0.20 and 0.25. The symmetry line is indicated in 
Fig. 4(c). The inset in (a) is the enlarged figure of the neigh- 
borhood of the Fermi level. 
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Fig. 4. The Fermi surfaces at n = 1.35 for (a) t' /t = 0.00, 0.10, 
0.15, (b) t'/t = 0.21, (c) t'/t = 0.22 and (d) t' /t = 0.23. 



level is located near the van Hove singularities, the shape 
of the Fermi surface is also sensitive to t' . With increasing 
t', the shape changes from Fig. 4(a) to Fig. 4(d). In Fig. 5, 
we display the bare spin susceptibility Xo{q, 0) along the 
symmetry line and the typical structure at t'/t = 0.23 
corresponding to Fig. 4(d). The introduction of t' leads to 
the remarkable wave-vector dependence: the incommen- 
surate peak or the broad peak structure near F point, 
and the hump structure near K point. The enhancement 
of xo(0, 0) at F point corresponds to the increase of the 
density of states by the van Hove singularity. On the 
other hand, the hump structure near K point originates 
from the peculiarity of the crystal structure. The wave- 
vector at K point in xo{q, 0) links not only F point and 
K point in the energy dispersion, but also the neighbor- 
ing K points. Thus, the hump structure also reflects the 
effect of the van Hove points around K point in . The 
peak structure at F point in Xo{q,0) implies the ferro- 
magnetic instability near t'/t ~ 0.21. Thus, the triplet 
superconductivity with Bi symmetry may be located in 
the neighborhood of the ferromagnetic phase, when we 
take the magnetic instability into account. 

Fig.6 shows A(fc) and J^{k) at n = 1.35, t'/t = 0.23, 




Fig. 5. (a) The bare spin susceptibility Xo(9jO) along the sym- 
metry line at n = 1.35 and T/t = 0.01 for t'/t = 0.00, 0.10, 0.15, 
0.20, 0.21, 0.22, 0.23, 0.24 and 0.25. (b) Xo{<l,0) at t'/t = 0.23 
and (c) the contour plot. We can see peak or hump structures 
around F and K points. 




Fig. 6. (a) The gap function A{k) with the /-wave triplet pairing 
(Bi) obtained at n = 1.35, t'/t = 0.23, U/t = 1.8 and T/t = 
0.01. (b) The contour plot of J^(fe) = \g{k)\^A{k). The absolute 
value is very large near K point on the Fermi surface. 



U/t =1.8 and T/t = 0.01. The symmetry of the gap 
function clearly possesses the /-wave {Bi) symmetry. 
The line node crosses each Fermi pocket around M point 
in Fig. 4(d). The absolute value of J^{k) becomes large 
near the Fermi surface around K point. This is because 
the Fermi velocity becomes small near the van Hove 
points. Thus, in this symmetry, the large superconduct- 
ing gap is open in the vicinity of the van Hove points, 
and the energy gain is large. Therefore, this solution is 
very stable. 

In the case of t' = 0, the eigen values are ~ 0.05 at 

most, and we can not expect superconductivity with high 
Tc- Rather, it means that the effect of the frustration in 
the ideal triangular lattice suppresses any instabilities. 
The introduction of t' is crucial to avoid such frustration. 
In Fig. 7, we display T-dependence of eigen values for 
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Fig. 7. The temperature dependence of eigen values for the /- 
wave pairing (Bi) at n = 1.35 and t' /t = 0.23; the thick line with 
black triangles for the third-order calculation, the thin line with 
white triangles for the second-order calculation and the thick line 
with black squares for the FLEX approximation. 



the /-wave triplet pairing at n = 1.35, t'/t = 0.23 and 
U/t = 1.8. For comparison, we have also evaluated the 
gap equation in the second-order perturbation expansion 
(Eqs.(3a) and (4a)), and in the FLEX approximation 
{U/t — 2).^'^'^'') In the both calculations, the same /- 
wave symmetry as discussed above is the most stable. 
The eigen values are illustrated together in Fig.7. It is 
interesting that the eigen value in the second-order cal- 
culation is about the same as that in FLEX. The pairing 
interaction in FLEX is given by the random phase ap- 
proximation of Xaio) composed of the dressed Green's 
function. In this case of U/t = 2, the Storner-type en- 
hancement factor 1/(1 — Uxo{<l)) is about 30. Simply 
speaking, this provides the pairing interaction 30 times 
larger than that in the second-order calculation. How- 
ever, this enhancement also yields a strong depairing 
effect by the normal self-energy. The strong depairing 
effect suppresses Tc too much. On the other hand, the 
third-order calculation provides finite Tc for the moder- 
ate U/t — 1.8. The third-order scattering process, es- 
pecially including the particle-paritcle ladder, weakens 
the depairing effect, which is over-estimated in FLEX. 
Simaltaniously, it provides the strong attractive interac- 
tion for the spin-triplet pairing. Thus, vertex correction 
terms other than simple one-boson propagation terms 
are important for the spin-triplet superconductivity, like 
the case in Sr2Ru04.^''^ In this case, the convergence 
of solutions with respect to U is very good. Both the 
second- and third-order processes cooperatively stabilize 
this solution. The transition temperature takes a reason- 
able value Tc ~ 0.005i ~ 5K for lOt ~ leV.^^^ Thus, this 
/-wave triplet pairing (i?i) is a possible candidate for the 
C0O2 superconductor. In fact, since Bi pairing state is 
a non-degenerate one with line nodes, different from 2D 
representations Ei and i?2, then we can naturally explain 
the line- node like behavior observed in 1/Ti. 

In conclusion, we have investigated possibility of un- 
conventional superconductivity in the repulsive Hubbard 



model on the 2D triangular lattice with the next neighbor 
hopping t'. We have evaluated the gap equation using the 
third-order perturbation expansion with respect to the 
on-site repulsion U. We have found that the /-wave spin- 
triplet pairing state {Bi) is the most stable in the wide 
range of t' and the electron density n more than half fill- 
ing. In this case, the introduction of t' is crucial to avoid 
the frustration in the triangular lattice, and move the 
van Hove points near the Fermi surface around K point. 
The bare spin susceptibility shows the broad peak around 
r point. These conditions stabilize the /-wave triplet 
pairing. From comparison with the second-order calcula- 
tion and the FLEX approximation, the third-order vertex 
corrections are important for the spin-triplet supercon- 
ductivity, rather than the remarkable ferromagnetic spin 
fluctuation. This is also the case in Sr2Ru04.") The /- 
wave triplet pairing (Bi) is the primary candidate for 
the superconductivity in Nao.35Co02.1.3H20, as far as 
the van Hove points near the Fermi surface arc located 
around K points. We can expect the case from the recent 
LDA band calculation.^^' In addition, the large relax- 
ation rate and uniform susceptibility suggest high den- 
sity of states at the Fermi level.'*' These support our 
theoretical proposal. At the present, in order to obtain 
a more reliable result about the possibility, we proceed 
with a further investigation using the band structure on 
the basis of the LDA band calculation. 
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